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We investigate the one-dimensional p- wave superconducting model with periodically modulated 
hopping and show that under time-reversal symmetry, the number of the Majorana zero modes 
(MZMs) strongly depends on the modulation period. If the modulation period is odd, there can be 
at most one MZM. However if the period is even, the number of the MZMs can be zero, one and two. 

In addition, the MZMs will disappear as the chemical potential varies. We derive the condition for 
the existence of the MZMs and show that the topological properties in this model are dramatically 
different from the one with periodically modulated potential. 
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Introduction. —Recently, searching for Majorana 
fermions (MFs) in condensed matter systems has at¬ 
tracted much attention,due to the following reasons. 
Firstly, MFs are their own antiparticles which were 
hypothesized to exist, but have long been elusive from 
detection. Secondly, because of their nonlocality and 
non-Abelian statistics, the zero-energy MFs, also called 
Majorana zero modes (MZMs), are proposed to be 
possible to realize fault tolerant topological quantum 
computation.® 7 Thus from both the fundamental and 
practical point of view, finding MFs is of great impor¬ 
tance. In condensed matter systems, MFs can appear as 
quasiparticle excitations in topological superconductors 
and MZMs refer to the zero-energy in-gap excitations. 
There are several suggestions of physical systems that 
may support the MZMs,®^ among which the one¬ 
dimensional p -wave superconducting (SC) model (also 
called the Kitaev model) P due to its simplicity and 
elegance, is the most studied one. Possible realization of 
the Kitaev model includes quantum wires with a strong 
spin-orbit coupling (or topologically insulating wires 
subject to a Zeeman mag netic field) and in proximity 
to a superconductor. 1112 J In addition, it can also be 
realized in cold-atom systems P 14 

Up to now, most of the theoretical works focus 
on ideal homogeneous 5 or potential-modulated Kitaev 
chainsf 15 17 or K itaev chains with longer-range hopping 
and pairingp ^ or even quasi-one-dimensional Kitaev 
chains with a finite width. 19 20 Particularly in the pe¬ 
riodically potential-modulated casep^HED it was found 
that under time-reversal symmetry, the number of the 
MZMs can be at most one and if the potential vanishes 
at certain sites, then the MZM will be very robust and 
stable for arbitrary strength of the modulation. How¬ 
ever, a very important problem unaddressed is the sta¬ 
bility and fate of the MZMs under hopping modulation. 
In this work, we investigate the hopping-modulated one¬ 
dimensional p-wave SC model which is an extension of 
the original Kitaev model. We found that, under time- 
reversal symmetry, the number of the MZMs strongly 


depends on the period of the modulation. If the period 
is odd, there can be at most one MZM. However if the 
period is even, in some parameter regimes the number of 
the MZMs can be two. Furthermore, the MZMs will dis¬ 
appear as the chemical potential varies no matter the pe¬ 
riod is odd or even. Therefore the topological properties 
of the hopping-modulated model are drastically different 
from those of the potential-modulated one. 

Method .—We consider a one-dimensional Kitaev p- 
wave SC model where the hopping is periodically modu¬ 
lated, the Hamiltonian can be written as 

H = J2(-ti4c i+1 + AcjCj+i + H.c.) + Vcjcj 

= \cm BdG C = Sj ] QQ ] H BdG QQ^C = I 

L 

= Y E NnVn, (1) 

n= 1 

where U = cos(27 via + S) is the periodically modulated 
hopping integral. A ^ 0 is the p-wave SC pairing gap 
and V is the chemical potential. Here a = p/q is a ratio¬ 
nal number with p and q being coprime integers. C ^ = 
(cI,...,4,Ci,...,c l ) and = (?7i,..., r] L , j?|,..., rf L ), 
with L being the number of the lattice sites. In addition, 

(~E l \ 

A = 

E 1 

V El) 

and Q is a unitary matrix that diagonalizes HbcIG- 
Defining Majorana operators and 'yf as 

if = 4 + c i)7 f = i(ci - 4) 
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then the quasiparticle operator r can be expressed as 

4 = ( 4 ) 

Z 2=1 

with (\) n ^ and ip n ^ being the amplitudes of the MFs Ti 4 
and in the nth eigenstate, respectively. If there exist 
MZMs, then none of the E n in Eq. 0 is zero under 
periodic boundary condition (PBC) while some of them 
become zero under open boundary condition (OBC) and 
the number of the MZMs is the number of the zero E n . 

Since U is modulated with a period q (the unit cell is 
enlarged by q times), therefore under PBC, we have 


L/q q -1 

H = + Ac s,Es+l,l) 

1=1 s =1 

—tq c \j c ld+l + ACqjCiJ+i + H.C. 

L/q q 

+ZZ w W ( 5 ) 

1 = 1 S=1 


Using Fourier transform, cj { = \/q/L^2 k c\ k e lkql , k G 
(—7r/g,7r/g]. In momentum space, we get 


H = 

*1 = 

H k = 
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M k \ 

At 


c q, — k\ 
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with the nonzero matrix elements of M k and A k being 
Mf ,s+1 = Mf +1,s = and A*’ s+1 = -A* +1 ’ s = -A, 
for s = l,...,g — 1. M k s = U for 5 = 1 , ...,g, 
while M^’ 1 = = —t q e~ lkq and A^’ 1 = —= 

—Ae~ ikq . 

Since we assume U and A in Eq. 0 is real (up to a 
global phase) throughout the paper, thus Hk respects the 
time-reversal, particle-hole and chiral symmetries and it 
can be unitarily transformed to an off-diagonal matrix 
apHU 

UH k U' = (// A 0 k ^j ,A k =M k + A fc , 

u = e ""' /4 = 75 (' '/) • < 7 > 


here is a Pauli matrix acting on the particle-hole space. 
Then the system belongs to the class BDI which is char¬ 
acterized by the Z index while the number of the MZMs 
can be represented by W which is calculated through 


w = =* r / 9 <H, 

7T Jk =0 Z k 

z _ e iOk _ D et (Ak) 
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FIG. 1: The energy spectra for o = 1/3 under OBC. Here 
A = 1 and L = 1632. (a), (b), (c) and (d) are for V = 0, 0.1, 
0.2 and 0.3, respectively. Only \E n \ < 0.05 are plotted. 



FIG. 2: (Color online) Det(A k= o)Det(A k=7T / 3 ) as a function 
of S and V, for a — 1/3 under PBC. The black solid, red 
dashed, blue dotted and green dash dotted lines are for V = 0, 
0.1, 0.2 and 0.3, respectively. The gray dotted line denotes 
Det(Ak=o)Det(A k=7V / 3 ) — 0. Here A = 1 and L — 1632. 


In fact, W just counts how many times the determinant 
of A k crosses the imaginary axis as k evolves from 0 to 
n/q. 

Results .—First we consider the a = 1/2 case, where 
t\ = —£2 = — cos S. Under PBC, we have 


M k 


A k 

Det(A k ) 


( V ~(ti+t 2 e 2ik )\ 

V-( t!+t 2 e- 2ik ) V )' 

( 0 -A(l-e 2ik )\ 

\A(1 - e~ 2ik ) 0 ) 5 

V 2 — t\ — t\ + 2A 2 — 2(£i£ 2 + A 2 ) cos 2k 
+2iA(£i T t 2 ) sin 2k 

V 2 + 4(A 2 — cos 2 S) sin 2 k. (9) 


Since Det(A k ) is real, W must be zero otherwise Det(A k ) 
will be zero for some k (which means the bulk energy gap 
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FIG. 3: (Color online) The distribution of the zero-mode 
MFs along the one-dimensional lattice for a = 1/3 under 
OBC. Here A = 1, L = 1632, V = 0.3 and 5 = 0.32 tt, as 
denoted by the red arrow in Fig. 0 d )- 
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FIG. 5: (Color online) The distribution of the zero-mode 
MFs along the one-dimensional lattice for a — 1/4 under 
OBC. Here A = 1, L = 1632, H = 0 and 8 = 0.253tt. (a) 
The eigenstate corresponding to E\ =0. (b) The eigenstate 
corresponding to E 2 — 0. 
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vanishes), therefore, there can be no MZMs. 

Generally, the periodic modulation can take many 
forms. For arbitrary G and £2 (G 7^ —G), we have 

Det(A k= 0) = V 2 - (h + t 2 ) 2 , 

Det(A k=v/2 ) = V 2 -(h-h) 2 +4A 2 . (10) 

If V 7^ 0 and Det(Ak=o)Det(A k=7r / 2 ) < 0? then Det(A k ) 
will cross the imaginary axis once as k changes from 0 to 
7 t/ 2. In this case, W m 1 and one MZM exists. Interest¬ 
ingly, at V = 0, we have 



Aik 

A^k 


( 0 A lk \ 

\A 2k 0 ) ’ 

— [ti + t 2 e 2ik + A(1 — e 2ik )}, 
-[t!+t 2 e- 2ik - A(l-e~ 2ik )}. 


( 11 ) 


In this case, the system can be divided into two separated 
subsystems A\ k and A 2k . As k evolves from 0 to 7r/2, we 
have 


Aik=oA lk=7T / 2 — (A + G) 2 — (A — G) 2 , 
A 2 k=oA 2k=7r / 2 = (A — G) 2 — (A + G) 2 - (12) 


FIG. 4: The energy spectra for a = 1/4 under OBC. Here 
A = 1 and L = 1632. (a) V = 0. (b) V = 0.1. Only 

\E n \ < 0.05 are plotted. 


Therefore, if both A lk =oA lk=7T /2 and A 2k =oA 2k=7r /2 are 
less than zero, then two MZMs will show up in this case 
and the existence of these two MZMs has been numer¬ 
ically verified (for example, G = 0-3, £2 = —0.8 and 
A = 0.5). 
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For a — 1/3, we have 



( V 


—tze 3lk 
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V 



• 

y— t^e 3tk 

-h 

V 


( 0 
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A 
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A 
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where t\ — cos(27r/3 + 5), £ 2 = cos( 47 r /3 + 8) and £3 = 
cos (5. In this case, 


Det(Ak) = [V (—3 + 6 A 2 + 2V 2 ) — cos 3 k cos 38 

+£A(4A 2 — 3) sin3fc]/2, (14) 


and Det(A k=0:7r/3 ) = [=fcos3(5 + F(-3 + 6 A 2 + 2 F 2 )]/2 
(— and + are for k = 0 and 7 r/ 3 , respectively). If 
Det(A k= o)Det(A k=7r / 3 ) < 0, then Det(A k ) will cross the 
imaginary axis exactly once as k evolves from 0 to 7 r/ 3 , 
which means that W = 1 and one MZM exists in this 
case [Ei in Eq. © is zero under OBC]. On the contrary, 
Det(A k=0 )Det(A k=7r/3 ) = 0 means that the bulk energy 
gap vanishes and Det(A k= o)Det(A k=7r / 3 ) > 0 means that 
Det(A k ) will not cross the imaginary axis as k evolves 
from 0 to 7 r/ 3 , in both cases W = 0 and no MZMs exist. 
Specifically, for 8 = (2ra + l)7r/6 with m being an integer, 
we have cos 3(5 = 0, therefore Det(A k=3 )Det(A k=7r / 3 ) ^ 
0 and no MZMs can exist, irrespective of the values of A 
and V. For example, we set A = 1 and L = 1632. In 
Figs, [l] and [ 2 ] we plot the energy spectra under OBC and 
Det(A k=3 )Det(A k=7r / 3 ) under PBC, respectively. As we 
can see, at F = 0, MZM exists for any 8 , except for 
8 = 7 r/ 6 , 7 r/ 2 , 57 t/ 6 , ?7 r/ 6 , 37 r/ 2 , 117t/ 6 where the bulk en¬ 
ergy gap closes. As V increases, for some 8, MZM van¬ 
ishes and as 8 evolves from 0 to 27r, topologically trivial 
(without MZM) and nontrivial (with one MZM) phases 
appear in turn. A typical distribution of the zero-mode 
MFs is shown in Fig. [3] and we can see that the two 
MFs and jf are well separated in real space and are 
located at the left and right ends, respectively while the 
actual decay length of these two MFs increases/decreases 
as the bulk energy gap decreases/increases. Finally when 
V > 0.32 where the condition \V(—3 + 6 A 2 + 2V 2 )\ > 1 
is satisfied, MZM disappears for any 8 and there is only 
topologically trivial phase and indeed for 8 = ( 2 m+l) 7 r /6 
with m being an integer, MZMs do not exist for any V. 

For a = 1/4, we have t\ = —t 3 = — sin (5 and £2 = 


—£4 = — cos (5. In this case, 


( V 

~tl 

0 -Ue Aik \ 

-ti 

V 

—£2 0 


0 

-t 2 

V -t 3 


t4e~ 4tk 

0 

-h v 

J 

( 0 

-A 

0 Ae 4ik \ 

A 

0 

-A 0 


0 

A 

0 -A 

5 

\-Ae~ Aik 

0 

A 0 J 



Det(A k ) = [3 — 8 A 2 + 8 A 4 + 8 ( 2 A 2 — 1)V 2 + 4V 4 
+ cos 4(5 + cos4fc (—1 + 8 A 2 — 8 A 4 
+ cos4(5)]/4. (15) 

At first glance, since Det(A k ) is real, thus, similar to the 
a = 1/2 case, there should be no MZMs. However, this 
is not true at V = 0. In the following we set A = 1 
and L = 1632 as an example. As we can see from Fig. 
[4j indeed at V = 0.1, there are no MZMs. However at 
V » 0, MZMs exist for 7 t/4 < 8 < 37t/ 4 and 57 t/4 < 
8 < 77 r/4. These MZMs are doubly degenerate [both 
Ei and E 2 in Eq. © are zero under OBC] and the 
distribution of the zero-mode MFs is shown in Fig. [5j 
The existence of these two MZMs can be explained as 
follows. At V = 0, we found that the eigenvalues of H k 
in Eq. © are doubly degenerate, therefore H k can be 
divided into two independent subsystems by a unitary 
transformation as 

PH “ P ' = sj ■ (16) 


here both Hi k and H 2k are 4x4 matrices while their 
eigenvalues are exactly the same. The unitary matrix P 
can be written as 


/I 

0 

1 



0 

-1 

\0 


1 0 0 

0 1 1 

-10 0 
0 1 -1 

1 0 0 

0 1 1 

1 0 0 

0 -1 1 


1 -1 0 

0 0 1 

1 1 0 

0 0 1 

-1 1 0 

0 0-1 
1 1 0 

0 0 1 


°\ 

-1 

0 

1 

0 

1 

0 

1 ) 


, ( 17 ) 


and 


Hj k — 
Det(Ai k ) = 


Mj k 
— A j k 


A j k A 

—Mjk) 


,4 = 1,2, 


[Det{A 2k )Y 

(A 2 - sin 2 (5) - (A 2 - cos 2 8)e 4ik . (18) 


At k = 0, Det(Ai k ) = Det(A 2k ) = cos 28 while at k m 
7t/4, Det(Ai k ) = Det(A 2k ) = 2 A 2 — 1 . If cos2(5(2A 2 — 
1) < 0, then both Det(Ai k ) and Det(A 2k ) will cross the 
imaginary axis exactly once as k evolves from 0 to 7 t/ 4, 
indicating that there exists one MZM in each subsystem 
and the number of the MZMs for the whole system is 
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two. Therefore for a = 1/4, at V = 0, the number of 
the MZMs are either two or zero while at V ^ 0, there 
are no MZMs. For general ti (i = 1,... ,4), at V / 0, 
Det(Ak) may not be real and there may exist one MZM. 
However at V = 0, the system can still be divided into 
two subsystems. In this case, if the conditions [(A + 
t 1 )(A + t 3 )} 1 2 3 4 5 6 7 8 9 10 11 12 - [(A -i 2 )(A-i 4 )] 2 < 0 and [(A-ii)(A- 
t 3 )] 2 — [(A + Z^KA+^j ] 2 < 0 are satisfied simultaneously, 
there will be two MZMs. 

Furthermore we found that, for general periodic mod¬ 
ulation, if the period q is odd, then the number of the 
MZMs is either zero or one. On the other hand, if q is 
even, then at F ^ 0, the number of the MZMs is still 
zero or one. However at V = 0, the system can always 
be divided into two independent subsystems and if the 
conditions 


[(A + t i)(A + £ 3 )... (A + t q - 1)] 2 
-[(A - t 2 )(A — U)... (A — t q )} 2 < 0, (19) 

and 

[(A — £i)(A — £3) ... (A — t q - 1)] 2 
— [(A + £ 2 )(A + £4).. . (A + t q )] 2 < 0, (20) 

are simultaneously satisfied, there will be two MZMs. 
However in this case, the subsystems cannot be further 
divided, making the maximal number of the MZMs be 
two in the whole system. 

Summary .—In summary, we have studied the num¬ 
ber of the MZMs and their stability in the hopping- 
modulated one-dimensional p-wave SC model. We found 
that the former strongly depends on the period of the 
modulation. If the period q is odd, there can be at most 
one MZM in the system while for an even q , the number 


of the MZMs can be zero, one and two. The existence of 
two MZMs can occur only at V = 0, since in this case, Ak 
in Eq. 0 can always be divided into two independent 
sub-matrices by a unitary transformation as 


- {L t) • 

s (j+t)/2,j = 1) for j = 1,3,5, 1, 

S U+Q)/‘2,j = 1, for j = 2,4,6, (21) 


At certain conditions, there exists one MZM in each sub¬ 
system and the number of the MZMs for the whole sys¬ 
tem is two. For the specific modulation form we consid¬ 
ered [ti = cos(27t ia + 8) with a = p/q\, only at q = 4n 
(n = 1,2,3,...) can Eqs. (19) and ( |20| ) be simultane¬ 
ously satisfied, therefore only in this case can there exist 
two MZMs. Furthermore, the MZMs will vanish as the 
chemical potential V varies. In the periodically potential- 
modulated model considered in Refs. nsmgandini when 
the time-reversal symmetry is present, there can be at 
most one MZM and if the potential vanishes at certain 
sites, then the MZM will be very robust and stable for 
arbitrary strength of the modulation. Clearly this is not 
the case in the periodically hopping-modulated model, 
therefore the topological properties differ drastically be¬ 
tween these two models. 
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